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ABSTRACT 

We examine the effect of a rapidly migrating protoplanet on a ring of planetesimals. 
Tlie eccentricities of the planetesimals are usually increased by Ae € (0.01, 0.1), with 
the exact increase being proportional to the protoplanet's mass, and inversely propor- 
tional to its migration rate. The eccentricity distribution is also substantially changed 
from a Rayleigh distribution. We discuss the possible implications for further planet 
formation, and suggest that the rapid passage of a protoplanet may not prevent the 
planetesimal disc from forming further planets. 
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1 INTRODUCTION 

A planet in orbit around a star will disturb the orbits 
of test particles. This situation, known as the 'restricted 
thre e body problem' h as been studied for centuries (see, 
e.g. iMurrav and Dermo tt (1999)). Particles close to reso- 
nant locations have their orbital elements changed due to 
repeated interactions with the planet. These changes may 
be computed using various analytic techniques (ibid), and 
the analysis can be extended to slowly migrating planets. 
Recently, a rapid migration m ode for prot oplanets has been 
discovered jMasset and PapaloizorJ I2OO3I) . The planet can 
have its semi-major axis halved in only a hundred orbits or 
so, which is far too fast for the effect on planetesimal orbits 
to be calculated analytically. 

In this paper, we attempt to quantify the expected ec- 
centricity increase in a planetesimal ring, due to a rapidly 
migrating protoplanet. A simple theoretical calculation is 
described in section |21 We describe a numerical model in 
section 13 and the results obtained from it in section 2] The 
possible implications for the formation of the Solar System 
are discussed in section |S] We summarise our findings in 
section |^ 



2 THEORETICAL PREDICTIONS OF 
PLANETESIMAL EXCITATION 

The orbits of the planetesimals can be be changed in two 
ways: 

(i) Long term interactions with resonances (mean motion, 
secular and co-rotation) 
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(ii) Close encounters with the planet 

Interactions between test particles and resonances have been 
studied for several centuries. However, for a rapidly migrat- 
ing planet, these resonances are much less important, since 
they will be sweeping through the disc with the migrat- 
ing planet. Hence, the resonance is unlikely to remain in 
the vicinity of a particle long enough to make significant 
changes. This is supported by the work of iTanaka and Idal 
(l99§), who found that a slow migration rate led to 'shep- 
herding,' (the planet gently brings the planetesimals with 
it), while fast migration leads to 'predatory' behaviour 
(the planet simpl y ploughs t hrough the particle disc) - see 
also the work of IWard and Hahn ( 1 9951) . In our calcula- 
ti ons, migration is even faster than the 'rapid' migration 
of iTanaka and Idal Jl999l) . 

iHaseeawa and Nakazawal il990l) found that the ex- 
pected change in eccentricity of a single particle per collision 
was given by 



(1) 



where Ri = 0.747, h = (m^i^^/Sm,)^^^ is the Hill parame- 
ter, and b is the impact parameter of the collision, in units 
of the semi-major axis. This is equation 38 of their pa- 
per, rewritten with conventional (rather than normalised) 
eccentricities. It is derived on the assumption that the stel- 
lar mass (m,) dominates, that the eccentricities and incli- 
nations are small, and that the impact parameter is large 
enough to ensu re that the particles can' t enter their mu- 
tual Hill sphere. iHasegawa and Nakazawal give another, sim- 
ilar, formula for the expected increase in inclination, but the 
coefficient is much smaller and we shall neglect inclination 
here (note also that the eccentricity increase is independent 
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of the inclination). Equation (with a different coefficient) 
may also be derived following the impulse approximation of 
iLin and Papaloizo Note the strong dependence on 

the impact parameter (curiously reminiscent of Rutherford 
scattering, although the link is not direct). This suggests 
that it is the closest encounter between the planet and plan- 
etesimal which will be the most important. 

But what is the closest encounter we should expect? 
Ever closer encounters will give larger changes (although 
note that equation eventually ceases to be valid), but en- 
counters with small impact parameters are rare. Both the 
planet and planetesimals are on near-circular orbits, so an 
impact parameter of zero implies that the two bodies are on 
the same orbit. Such particles would have a (near) infinite 
synodic period. Consider a particle at semi-major axis ao, 
and the planet migrating from ao -I- Aa to ao — Aa. We can 
expect an interaction with b ~ Aa/ao if the time for the 
planet to migrate between the two limits (2Aa/d) is equal 
to the synodic period for orbits at oq and ao-|- Aa. Formally, 
this neglects the orbital inclination, but we expect the effect 
to be small in a disc. We have 



rayn(ao, Aa) 
2a 



cross 

2Aa 



3Aa 
Aa 

ao 



kcp 



3ao 



2kcp(ao) 



(2) 
(3) 
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where v^^^^ is the circular Keplerian velocity, and Ty.^^ is the 
Keplerian orbital period. Substituting into equation ^ we 
find 



(Ae^)5 = 



Tkep(«0) 



ao 



(5) 



For a 50 M0 planet in orbit around the sun and migrat- 
ing at 10* cms" through a ring close to 3.5 au,^ equation |K| 
suggests (Ae^) 2 r; 0.03. Several approximations were made 
in deriving equation |K| It assumes a single close encounter, 
and the estimated timescale for this (based on the synodic 
period) is rather imprecise. The numerical factor could eas- 
ily be incorrect. At this point, numerical simulations become 
useful in determining the evolution. 



3 MODEL 

To simplify the system as much as possible, we considered 
the interaction of a migrating planet with a narrow ring of 
test particles. 

The bodies in our simulation belonged to one of three 
types: 



(i) The central star, m» 

(ii) The migrating planet, rripi^^ 

(iii) The test particles (planetesimals). 



''part 



We evolved our system using Newtonian gravity, plus a 
torque to migrate the planet. To avoid the computational 
expense of a full n-body calculation, the test particles did 



not interact with each other. However, the gravitational 
interactions between the star and planet, the star and the 
particles, and the planet and the particles, were all fully 
computed. 

We neglected the effect of gas damping, since we were 
interested in the effects of rapid migration. This is reason- 
able, since our largest migration timescales were still far 
shorter t han the d amping timescales. Using the notation 
of Tanak a and Idal (il99a) . we typically had Tmig < 1 while 
T'gas ~ 1000, even with enhanced gas density. This holds even 
if the gas damping timescale is assumed to be due to wave 
excita tion, as will be the case for the more massive planetes- 
imals iArtvmowic dllQQ^l . Hence we did not expect the gas 
damping to be sufficient to make our migrating protoplanets 
'shepherd' the planetesimals. 

3.1 Orbital Migration 

Orbital migration of the planet is imposed by applying a 
torque to the orbit of the star and planet. In the interests 
of simplicity, we apply a torque chosen to give a constant 
migration rate, d. The required torque may be computed 
using Kepler's Laws: 
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(6) 



where L = /la^Q, is the orbital angular momentum, and G 
is the required torque. The torque was applied as an extra 
force to the motion of the star and planet. Migration was 
halted once the planet's semi-major axis dropped below a 
preset value. 

3.2 Intial Conditions and Integration 

The star (IMq) and planet were positioned in a circular 
orbit about their centre of mass. The initial orbital sepa- 
ration was usually 6 au, and migration stopped at 0.5 an. 
We distributed the test particles in circular orbits around 
the centre of mass of the star-planet system. The particles 
were uniformly distributed in azimuth and radius (usually 
3 — 4 au). Finally, we added small perturbations to their cir- 
cular motion, following the prescription of Ste wart and Ida| 
i2000l) . to give a Rayleigh distribution of e and i (ibid). We 
tested a variety of planetary masses and migration rates. 

We integrated the equations of motion using a Runge- 
Kutta integrator with adaptive step-size. To avoid catas- 
trophically small timesteps, test particles which strayed too 
close to the star or planet were eliminated. Similarly, parti- 
cles which reached large radii were removed. At the end of 
each timestep, we computed the values of e, i and a for each 
particle using the Laplace-Runge-Lenz formalism. 



4 RESULTS 

4.1 General Behaviour 

When the planet was far from the planetesimal ring, waves 
in (e, a) space were observed propagating through the ring. 



^ Note that 10* cms" ^ 0.2auyr"^, which is n ot unreason- ^ For the short timescales we are studying this is reasonable, since 

able fo r the rapid migration discussed bv lMasset and Papaloizou the close proximity of the (relatively) massive planet will have a 

^2003^ far greater effect 
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These were caused by resonances sweeping through the ring 
- if the planet did not migrate, then sharp peaks in the ec- 
centricity were observed, corresponding to resonant orbits. 
However, the increase in eccentricity prior to the planet's en- 
counter with the ring was fairly low. As the planet ploughed 
through the ring, the planetesimals were catapulted up 
lines of constant Jacobi energy, Ej , given by llHavashi et alJ 
119771): 



■ 2 , .2\ 

e +1 ) 



8 



(a' 



(7) 



where a' is the particle semi-major axis in units of the 
planet's semi-major axis (this energy is expressed in scaled 
units). The encounters with the planet were generally at a 
distance further than the L2 point (this is consistent with 
the assumptions made in deriving equation 0. 




Figure 1. Two sample particle trajectories in e-a space 



4.2 Numerical Results 

Several runs were performed, with the parameters and re- 
sults summarised in Table For all these runs, the ring 
of particles lay between 3 and 4 au initially. In this ta- 
ble — + i^, the conventional measure of deviation 
from circular, coplanar orbits. In practice, ^ was dominated 
by eccentricity, not inclination (cf Section We used a 
Kolmogorov-Smirnov test to investigate if the final distri- 
bution of ^ was still drawn from a Rayleigh distribution. In 
every case, this hypothesis was strongly rejected.'' Qualita- 
tively, we found that there was a tail of high ^ particles. 
This made the RMS values for ^ somewhat changeable, so 
we quote the median value as well. For the remainder of 
this discussion, references to the final ^ values refer to the 
median. 

Most of the results of Table bear out the qualitative 
behaviour suggested by equation|S] but differ quantitatively. 
The final median ^ values are close to those predicted (recall 
that ^ was dominated by e in our simulations). Increasing 
the migration timescale or the mass of the migrating planet 
increases the final value of However, this doesn't quite 
happen in the linear manner predicted by equation |S] The 
detailed values are similar to those predicted, but not iden- 
tical. 

Figure0shows part of the problem, plotting two sample 
particle trajectories in e-a space. These particles have moved 
around significantly in e-a space, and have not done this just 
once. The trajectories are characterised by large jumps (cor- 
responding to those analysed above) interspersed with pe- 
riods where the particles are almost stationary. Sometimes, 
the eccentricity is pumped up more slowly as well (presum- 
ably when the particle happened to be close to one of the 
resonances of the migrating planet). This is in sharp con- 
trast to the 'single encounter' approximation of Section |5| 
Obtaining better predictions of the final ^ value is therefore 
problematic. Even if the resonant pumping is neglected, we 
cannot apply a simple random walk approach. Firstly, the 
number of steps (large 'jumps') is also random - and not 
large. Worse, the steps aren't truely random, but are con- 
strained to have constant Jacobi energies (although the con- 

^ We also checked that the KS test did allow the initial distribu- 
tion to be Rayleigh 




Figure 2. Surface of median ^ values following migration 

stant value changes for each step, since the planet will have 
migrated in the meantime). 

Run 8 behaved slightly differently. However, even this is 
not surprising when compared to run 7. The two runs were 
identical apart from the initial ^ value, and the initial ^ for 
run 8 was higher than the final value for run 7. It is therefore 
not surprising that the final ^ value for run 8 is higher than 
that for run 7 - but notice that the increase in ^ is similar 
for both runs. 

We also ran a grid of 40 models, with rUpi^^ ranging be- 
tween 10 and 400 Mgj and a in the ra nge 2x10^ to lO'^ cms"^. 
This covers the range expected by iMasset and Papaloizoul 
(2003) to undergo runaway migration, and a bit more on 
each end. All these runs started with average e and i values 
of 10"^ Fi gure|5| plots a surface showing the resultant me- 
dian 5 values. Planets more massive than 100 Mq frequently 
managed to increase ^ to be greater than 0.1 (and reached 
0.4 for a 400 M^ planet migrating at 2 x 10^ cms"^). Al- 
though not a perfect match, the behaviour predicted by 
equation|S]is seen. 



5 IMPLICATIONS 

One major problem with the migration of planets in their 
nascent discs is t hat protoplanets "coul d be too mobile for 
their own good" JWard and Hahnll2000l) . The rapid migra- 
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Table 1. Numerical Results. The quantity f is defined as = 



Recall that 10* cms"! fa 0.02 auyr'l 



Run 


plan' a? 


a/cms ^ 


Initial e,i RMS 


Final ^ RMS 


Final ^ median 


1 


100 


10"* 


10-3 


0.0859 


0.0360 


2 


100 


5 X 103 


10-3 


0.1472 


0.0978 


3 


100 


10* 


2 X 10-3 


0.0931 


0.0358 


4 


100 


5 X 10^ 


2 X 10-3 


0.1717 


0.0987 


5 


200 


10-1 


10-3 


0.1549 


0.0839 


6 


10 


10* 


10-3 


0.0238 


0.0039 


7 


10 


5 X 10^ 


10-3 


0.0612 


0.0071 


8 


10 


5 X 10^ 


10-2 


0.0352 


0.0163 


9 


10 


2 X 10^ 


10-3 


0.0411 


0.0180 



tion mode of lMasset and Papaloizoul (|200^ makes the pro- 
toplanets even more mobile - it seems to be very easy to lose 
protoplanets into the Sun. Whilst this is obviously not very 
'good' for the protoplanets in question, would such events 
make the formation of the Solar System impossible? Put 
another way, "Is the current Solar System the first Solar 

System?" 

Considering the Q parameter of iToomrel il964l) sug- 
gests that the disc around the young Sun may have been 
up to five times more massive than the Minimum Mass 
Solar Nebula ( MMS N) model of Hayashi (1981) (although 
iHavashi et al.l il985ft is likely to be a more accessible refer- 
ence) . This would give plenty of material for forming several 
generations of planets. The question is therefore whether a 
migrating planet will disrupt the disc sufficiently to prevent 
further planet formation within the nebula lifetime of a few 
Myr. 

If the migration rate is low enough to permit shepherd- 
ing (and hence depletion o f planetesimals in ward of the pro- 
toplanet's intial location'). lArmitaed i2003l) has shown that 
it is unlikely that fresh material could diffuse into the inner 
portions of the disc. In this case, further planet formation 
would not be possible. We have considered a migration rate 
too high for shepherding, so we must examine whether our 
planetesimal disc is too hot to allow further planet forma- 
tion. 

In our simulations, Ae is typically in the range 
(0.01,0.1), dependent on the mass of the protoplanet and 
its migration rate. Planetesimals are fairly fragile objects - 



bodies with m 



msl- 



10 K will suffe r disruptive collisions 



if e > 0.01 jKobavashi and IdallioOlh . However, gas damp- 
ing can recircularise the orbits of such bo dies very rapidly 
- prob ably less than lO"* yr (see figure 1 of iTanaka and Ida 
il999t) and references therein) - especially if the disc is denser 
than the MMSN. The damping timescale will reach a peak 
for rripart ~ 10^^ g, but these bodies require e > 0.1 before 
sufi'ering disruptive collisions. Still larger bodies will require 
even higher e values before they shatter, but these damp 
faster due to the excitation of density waves dArtymowicd 
^93). Gravitational focussing will also fall as eccentricities 
increase, due to the higher relative velocities implied. This 
w ill cause a dramatic fa l l in th e collision rate (cf figure Al 
of lWeidenschilling et al.l (^23)), so the likelihood of disrup- 
tive collisions is reduced. Finally, in a real planetesimal disc 
there will be a distribution of sizes. Wcthcrill and Stc warti 
il993 f) found that grinding moderately sized planetesimals 
into rubble helps the larger planetesimals accrete them, so 
some excitation may even be helpful. 

Based on this discussion, it seems that the rapid passage 



of a protoplanet will not prevent further planet formation. 
Coagulation will be suppressed for a while, but this should 
be short compared to the disc lifetime. However, the ar- 
guments in the preceding paragraph are not rigorous, and a 
longer term calculation of the 'end' states would be required 
to give a firm answer. This is particularly true just as the 
planetesimal disc re-achieves equilibrium. 



6 CONCLUSIONS 

Combining our results (including some varying no not ex- 
plicitly listed here), we conclude that a rapidly migrating 
protoplanet will increase the eccentricity of a planetesimal 
ring by roughly 



7 X 10' 



''^plaii 



ap 
1 an 



(8) 



where ao is the initial semi-major axis of the planetesimal 
ring and a is the protoplanet's migration rate. During the 
passage through the ring, the distribution of eccentricities 
will also be substantially changed from a Rayleigh distribu- 
tion. In particular there will be a tail of high eccentricity 
bodies. 

Simple arguments suggest that this effect should not 
prevent subsequent planet formation by the surviving plan- 
etesimals (very few are ejected or accreted). The resultant 
eccentricities should not be dangerously high for at least 
some of the planetesimals, collisions are less likely due to re- 
duced gravitational focussing, and damping timescales are 
short. However, problems could arise shortly before the plan- 
etesimal disc achieves equilibrium once more, and further 
simulations (spanning longer timescales and including all 
gravitational interactions) are needed to assess the effect in 
detail. 
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